Rectification of thermal fluctuations in a chaotic cavity heat engine 
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We investigate the rectification of thermal fluctuations in a mesoscopic on-chip heat engine. The 
engine consists of a hot chaotic cavity capacitively coupled to a cold cavity which rectifies the 
excess noise and generates a directed current. The fluctuation-induced directed current depends 
on the energy asymmetry of the transmissions of the contacts of the cold cavity to the leads and 
is proportional to the temperature difference. We discuss the channel dependence of the maximal 
power output of the heat engine and its efficiency. 
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I. INTRODUCTION 

Rectification is central to the operation of electrical 
circuits. More than 60 years ago, Leon Brillouin, then 
at IBM, raised the question of whether an electric cir- 
cuit consisting of a resistor and a diode can become a 
Maxwell demon rectifying its own thermal fluctuations. 1 
Using inappropriate generalizations of Langevin dynam- 
ics for systems with nonlinear diffusion coefficients could 
indeed lead to such rectification, in obvious violation of 
the second law of thermodynamics. 2 Brillouin's paradox 
was solved by taking into account the diode's contact po- 
tentials. 3 Later, it was shown that a system with diode 
and resistor at two different temperatures cannot exceed 
Carnot efficiency 4 in agreement with the second law. 

Nowadays, thermoelectrics is of increasing importance. 
In the continuing quest for smaller scale electric circuits 
the evacuation of heat proves to be a major obstacle. 
Therefore it is interesting to explore whether some of the 
energy that is dissipated can be harvested and put to use. 

There are many ways to generate directed currents. 
In recent years Brownian particles in ratchets subject to 
periodic driving have found much interest in very differ- 
ent fields of science. 5 Here we are concerned with a more 
subtle form of driving: The only external agent acting on 
the system is noise that can be generated by an external 
thermal equilibrium bath. External noise can generate 
directed currents even in periodic potentials with inver- 
sion symmetry if the noise power depends on the location 
of the Brownian particle. 6 10 Such state dependent dif- 
fusion is also at the origin of the difficulties encountered 
by Brillouin. 1 

In this paper, we investigate the rectification of ther- 
mal fluctuations into a directed electric current in a meso- 
scopic heat engine. The latter consists of two capacitively 
coupled chaotic cavities arranged in a three-terminal ge- 
ometry as shown Fig. 1. The upper cavity is the rectifier 
and is connected via two contacts with energy-dependent 
transmissions to electron reservoirs. The lower cavity 
provides the external source of thermal noise. It is con- 
nected via only a single contact to another electron reser- 




FIG. 1. (Color online) Schematic of the double cavity and 
operation scheme: Electrons enter the cavity, gain an energy 
eSU and then leave the cavity again, experiencing different 
transmissions of the contacts. 



voir. The energy dependence of the contact transmissions 
is a generic feature of mesoscopic conductors and leads 
to an intrinsic nonlinearity in the upper cavity. 11 The 
three-terminal setup allows for separated heat and charge 
currents in contrast to two-terminal setups where these 
currents are necessarily aligned. 

The mechanism giving rise to the current is shown 
schematically in Fig. 1. Electrons enter the cold up- 
per cavity at an energy E. They absorb the energy eSU 
from the fluctuating potential generated by the hot lower 
cavity in the upper cavity and afterwards leave the cold 
cavity again. As the transmissions through the upper 
cavity's contacts are energy-dependent, the ratio between 
transport processes involving the left and right lead is dif- 
ferent at energies E and E + e5U, thus giving rise to a net 
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electrical current through the cavity. The rectification is 
controlled by Coulomb coupling which is dominant at low 
temperatures, whereas at higher temperatures phonon ef- 
fects become important. 12 

Thermoelectric properties of both open mesoscopic 
cavities 13,14 and in the Coulomb blockade limit 15 17 have 
been of interest. The quantization of energy levels in 
small dots leads to exceptional thermo-electrical prop- 
erties. 18 Both in two- and three-terminal structures the 
limit in which the ratio of electric to heat current is given 
by the ratio of the charge to an energy quantum can 
be reached. 19,20 Still, while the efficiency of such nano- 
engines can in principle be optimal leading to an infinite 
figure of merit ZT 21 , the current they deliver is small, 
typically of the order of 1 pA. It is therefore of interest 
to explore how power output and efficiency scale as dots 
are opened and turned into cavities with contacts that 
permit currents that are much larger than the tunneling 
current of a Coulomb-blockaded quantum dot. 

The physics of Coulomb coupled conductors is of inter- 
est in nanophysics for on-chip charge detectors, 22 quan- 
tum Hall edge states, 23 and the Coulomb drag in which 
one system that carries a current induces a current in a 
nearby unbiased conductor. 24 27 In the setup of Fig. 1 the 
current carrying conductor in the Coulomb drag problem 
is replaced by an unbiased but hot conductor. 

Our paper is organized as follows: In Sec. II, we present 
our model of the double cavity and describe our theoret- 
ical approach. Our results are presented in Sec. Ill and 
conclusions are drawn in Sec. IV. Calculational details 
are presented in the Appendices. 



II. MODEL AND METHOD 

We investigate transport through two capacitively cou- 
pled open quantum dots with mutual capacitance C, cf. 
Fig 1. Each cavity z = 1,2 is coupled via quantum point 
contacts (QPCs) to electronic reservoirs r = L, R. The 
latter ones are assumed to be in local equilibrium and 
described by a Fermi distribution fi r (x) = {exp[(x — 
Mir)/^B©i] + 1} -1 with temperature and chemical po- 
tential \±i r . Interaction effects are captured by capacitive 
couplings Ci r between the cavities and their respective 
reservoirs that leads to screening of the potential fluctu- 
ations. 

We consider the system in the semiclassical limit where 
the number Ni r of open transport channels 28 30 in the 
QPCs is large under conditions at which dephasing de- 
stroys phase information but preserves energy. We can, 
thus, characterize the chaotic cavities by a distribution 
function fi(E) that depends on energy only, and focus on 
a semi-classical description of the physics without coher- 
ence. 31 For later convenience, we write the distribution 
function as 

/i= y^ /ir +«y/i. a) 



Here, the first term describes the average value of fi and 
is given by the average of the distributions of the reser- 
voirs weighted with the transmission T ir of the respective 
QPC. The second term Sfi describes fluctuations of fi 
around its average. Additionally, each cavity is charac- 
terized by its potential Ui which also fluctuates by SUi. 

We assume the transmissions to be energy dependent 
which we model to first order as T ir = Tf r — eT' ir SUi. 
The energy dependent transmission leads to a nonlin- 
ear current voltage characteristic. Even without external 
noise such a nonlinearity requires a self-consistent treat- 
ment. 32,33 In our system, Fig. 1, we need a self-consistent 
treatment not only of the average Hartree potential but 
in addition the fluctuating potentials. We remark that 
while the energy-independent part Tf r scales linearly with 
the number of open transport channels A^ r , the energy- 
dependent part T[ r is independent of N ir . 

The starting point of our theoretical investigation is 
a kinetic equation for the distribution functions fi (see, 
e.g., Ref. 34), 



where z^p denotes the density of states of cavity i. It 
describes the change of charge in a given energy inter- 
val due to changes in the potential Ui, in- and outgoing 
electron currents through the QPCs as well as their fluc- 
tuations Six where the index E indicates summation over 
all contacts r of cavity i. 

Expressing the charge inside the cavities via the distri- 
bution functions fi as well as via the capacitances and po- 
tentials, we obtain a relation between Sfi and SUi which 
allows us to transform the kinetic equation (2) into a 
Langevin equation for SUi. Neglecting terms that are 
cubic and higher in the potential fluctuations, the latter 
can be converted into a nonlinear Fokker-Planck equation 
with a diffusion function that depends on the cavity po- 
tential. The nonlinearity of the Langevin equation leads 
to subtleties in the interpretation of the stochastic inte- 
gral (known as the ltd and Stratonovich problem) that 
gives rise to different Fokker-Planck equations. The "ki- 
netic prescription" of Klimontovich 35 provides a steady- 
state solution of Eq. (2) that is in global thermal equilib- 
rium, thus avoiding the Brillouin paradox mentioned in 
the introduction. 

We stress that the nonlinearity is the technical ori- 
gin of rectification in the cavity. A linear system would 
not exhibit this feature. From the Klimontovich- Fokker- 
Planck equation we obtain (SUi) and (SUiSUj), cf. Ap- 
pendixes A-D for details. The charge currents between 
upper cavity and contact r are given by 

hr = ^j dET lr (f lr - h) + SI r . (3) 
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FIG. 2. Inverse effective capacitance as a function of the 
capacitance between the two cavities C for different values of 
e 2 UF/C^. Note the different scale for the various curves. 



III. RESULTS 

The critical nonlinearity of the cavity is quantified 
by the amount of symmetry-breaking in the energy- 
derivatives of the transmissions of the upper cavity, given 
by the rectification parameter A: 



A 



G 1L G iR G 1R G 1L 



(4) 



where G ir = (e 2 /h)T° r and G' ir = (e 3 /h)Tl r . We will 
see below that it is the rectification conversion factor be- 
tween energy and charge for an unbiased cavity. This 
parameter appears in many places in terms of interac- 
tion corrections. For example, if we consider the un- 
coupled upper cavity (G = 0, i = 1), the single cav- 
ity conductance G\ = GilGir/Gi^ (which is the series 
combination of the left and right leads) has an interac- 
tion correction of AGi = — 2Gi /i G$]A 2 /c j BOi/G 2 £, where 
Cie = Y^ r C\r is the total capacitance of the upper cav- 
ity and G-^ 1 = (e 2 ^ip) -1 is its electrochemical ca- 
pacitance. We note that while Gi scales with the chan- 
nel number iV, the correction AGi is of order ex N~ x 
(quantum corrections in a coherent cavity 36 are of order 
oc N°). The rectification parameter A also appears in 
the second-order conductance. 

We now turn to the coupled cavities, cf. Fig. 1. Ap- 
plying different temperatures 0i and 02 to the reservoirs 
that couple to cavity 1 and 2, respectively, while keep- 
ing the electrical contacts grounded (Vi r = 0), we find to 
leading order in the nonlinearity a charge current through 
the cavity given by 



</lL> = — M©1 - ©2), 
TRC 



(5) 



where we assumed identical capacitances and densities 
of states for the two cavities. Now, we give a physical 
interpretation of each term in Eq. (5), the rectified cur- 
rent, trc — C e ff/G e ff denotes an effective RC time of the 



double cavity. It is determined by the effective conduc- 
tance of the double cavity, G e ff = Gi^G^s/ '(Gie + G^s) 
which is largest if both cavities have equal conductances. 
Furthermore, it depends on the effective capacitance 
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(6) 



describing how strong the interaction is between the two 
cavities. It should be minimized (without entering the 
Coulomb blockade regime) to maximize the rectified cur- 
rent. It grows as C~ 2 and for large couplings approaches 
the constant value (2G^ — G /i )G^/G /i , cf. Fig. 2. Next, 
as stated, the rectified current (5) is proportional to A, 
which characterizes the asymmetry of the system: The 
system is asymmetric if either the left-right conductances 
and/or their energy derivatives differ. Finally, the cur- 
rent (5) is linear in the applied temperature difference, so 
the rectified current is zero in global thermal equilibrium, 
as must be the case in order to satisfy the second law of 
thermodynamics. We note that the sign of the current 
flips under either exchange of the system lead nonlinear- 
ity or under exchange of the cavity temperatures. 

As the energy-dependent part of the transmission does 
not scale with the number of transport channels, the cur- 
rent Eq. (5) also turns out to be independent of the chan- 
nel number. For realistic values 28 30 of G e ff = 10 fF, G' = 
(e 2 //i)(mV) _1 and 6 2 - 6i = 1 K, we find / - 0.1 nA 
which can be readily detected in current experiments and 
is two orders of magnitude larger than currents through 
typical Coulomb-blockaded dots. 

In order to convert the heat extracted from the hot 
reservoir into useful work, we have to make the noise- 
induced current flow against a finite bias voltage Vll — 
Vir. The bias induces a counterflow of current given 
to leading order in the nonlinearity by Gi(Vll — Vir), 
thus reducing the total current. At the stopping voltage, 
Vstop = A /cb(Oi — @2)/(G\Trc), there is no current flow- 
ing through the system. The output power is given by 
P = (Iil)(Vil — Vir). It is parabolic as a function of the 
applied voltage, vanishing at zero bias and the stopping 
voltage; it has a maximum at half the stopping voltage 
given by 



P — 

J- max — 



A 2 



(fc B (0i-©2)) 2 



(7) 



Energy is transferred between the cavities in the form 
of dissipated power in the upper cavity, P = ^2 r (Ii r Vi r ) , 
the heat current given up by the lower hot cavity to the 
upper cold cavity, Jh = (^2^2)? an d the heat current 
given up by the upper cold cavity to its heat reservoirs, 
Jc = ^2 r ((Ui — V\ r )Ii r ). It is straightforward to check 
that Jh = Jc + P in our model, so energy is conserved 
in the system. The efficiency, 77, of the heat-to-charge- 
current converter is given by the ratio between the output 
power, P to the inter-cavity heat current, Jh- To leading 
order in the energy-dependent transmissions, this heat 
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current is given by 
Jh 



trc 



-fc B (e 2 -9i), 



(8) 



because heat will flow from hot to cold even without the 
nonlinearity. The correction to this result that is linear in 
the voltage applied across the upper cavity is suppressed 
by G' lr (as it must to satisfy an Onsager relation, see Ap- 
pendix E). As indicated earlier, the asymmetry parame- 
ter A controls the process of energy-to-charge conversion, 
(Iil)/Jh. 

The efficiency 77 = P/Jh exhibits the same parabolic 
bias dependence as the output power since the heat cur- 
rent is independent of the applied bias. Hence, for a given 
temperature difference, the maximal efficiency occurs at 
maximum power and is given by 



A 2 



Me 2 -ei). 



(9) 



For G\ = 5e 2 /h and parameters as above, we estimate 
P max ~ 2fW and a maximal efficiency of ~ 1% of the 
Carnot efficiency for a device working at liquid-helium 
temperatures. We note that while the maximal power 
of the system scales inversely with the number of avail- 
able transport channels, the maximal efficiency even de- 
creases with the number of channels squared. This is 
because for a large number of open channels the effect 
of the energy-dependence of the uppermost channel be- 
comes less important. This effect can be seen in Fig. 3 
where the logarithmic derivative of the QPC conductance 
is plotted, which controls the stopping voltage and other 
rectification figures of merit in the case where one con- 
tact is energy-independent. For a stronger nonlinearity, 
such as a truly step-like transmission, a nonperturbative 
analysis is required which could give rise to much higher 
efficiencies. 4 

In order to demonstrate that it is the nonlinearity of 
the rectifying cavity that is the key ingredient to our 
heat engine, we now briefly consider an alternative setup 
where a rectifier is coupled to a resistor R at temperature 
82 with a capacitor Cr in parallel instead of a second 
cavity (see Fig. 4 for a circuit diagram). Repeating the 
same analysis as above, we find that the heat-induced 
current, the maximal power as well as the efficiency at 
maximum power are given by the same expressions as 
above. The only difference is that the effective conduc- 
tance and capacitance take on different values, G e ff — » 
Gie, C eff [CrCv + C(C R + C E )][C M C E (C + cy + 
RG*{CC R C* + (C + C R )Cl - CC R C^}/(C 2 C^). 

We finally apply our results to the semiclassical regime 
in the limit G2S ^> Gis for simplicity. For large energies 
the conductance of the QPCs is given by 



G r (E) 



e 2 (E - E r 



h 



A? 



-G(E - E r 



(10) 



with E r being the energy that marks the transition from 
tunneling to ballistic transport and A r describing how 
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FIG. 3. Conductance of a QPC for both fully quantized 
steps (full black) and semiclassical interpolation (full red). 
Its logarithmic derivative (black and red dashed) controls the 
heat-to-charge conversion, stopping voltage, power and effi- 
ciency of the energy harvester. 



open the contact is. For equal conductances of the two 
QPCs at the Fermi energy, Gl(Ef) = Gr(Ef), we ob- 
tain for the maximal power 



e 4 7 2 [M@i - 6 2 )] 2 (1 - 1Z)\E F - E L ) 



7-2 



8h 



(ii) 

with 1Z = Al/Ar. We thus see that in order to maxi- 
mize the power, we need a strong asymmetry in the con- 
tacts, Al <C Ar while keeping the conductance of each 
contact the same. For the semiclassical result 7 = 1/2, 
we find that the power drops upon increasing the energy 
transport window E F — E^. We note, however, that for 
7 > 2 the contribution from the conductance will out- 
weigh the contribution from the nonlinearity and, thus, 
lead to a maximal power that increases with the trans- 
port window. While the efficiency will drop with the 
inverse square of the transport window, we remark that 
a conductance that is exponential in energy will have an 
efficiency that is independent of the window size. 



IV. CONCLUSIONS 

We have examined a mesoscopic energy harvester con- 
sisting of a pair of quantum dots and find that as the 
contacts are opened, the power output can increase but 
typically with a drop in efficiency for the weak nonlin- 
earity considered here. Our work demonstrates the im- 
portance of the asymmetric energy dependence of the 
contact transmissions. Energy harvesting from environ- 
mental fluctuations is an important goal. It might lead to 
nano-scale devices which can function independently of 
an external power supply. In densely packed electronic 
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FIG. 4. Circuit diagram of a cavity coupled to a hot resistor. 



circuits energy harvesting might alleviate the heat re- 
moval problem. Our results are useful for future experi- 
ments that realize solid state energy harvesters. 



ACKNOWLEDGMENTS 

We acknowledge support from the project 
NANOPOWER (FP7/2007-2013) under Grant No. 
256959. AN J acknowledges support from NSF Grant 
No. DMR-0844899 and the University of Geneva. RS 
was supported by the CSIC and FSE JAE-Doc program, 
the Spanish MAT20 11-24331, and the ITN Grant No. 
234970 (EU). 



Appendix A: Kinetic equation 

Our starting point is the set of kinetic equation for the 
distribution functions fa of the cavities, 



dfo dfi_ 
dt ~ dUi 



Ui 



7^E T *-/J + -^ ( A1 ) 



a constant part that is given by the average of the Fermi 
functions of the reservoirs weighted with the transmission 
of the respective QPC and a fluctuating part Sff. 



irjir 



y] ~FT~fir - AiLR(/iL - /»r) ( SUi 



Sfi 



(A2) 



Here, in the last step, we used T ir = Tf r — eT' ir SUi and 
expanded the whole expression up to second order in SUi. 
We, furthermore, introduced the asymmetry parameter 



A 



G' T GiR - G- R GiL 



iLR 



(A3) 



and abbreviated G ir = (e 2 /h)T? r , G' ir = {e^/h)T[ r , 
where r = L, R, £ refers to the left QPC, the right QPC 
or the sum over all QPCs adjacent to cavity i = 1,2. 

In order to relate the fluctuating part of the distribu- 
tion function Sfi to the fluctuating part of the potential 
SUi, we express the charge Qi C inside cavity i once in 
terms of its distribution function and once in terms of 
the potentials and capacitances, 

= eu iF J dE f^U + Sf^j - e 2 v iF U u (A4) 

Q*c = C *r( U i ~ Kr) + C ig (Ui ~ V ig ) + C(Ui ~ Ui)> 

(A5) 

where i denotes the index opposite to i. Equating the 
fluctuating parts of both equations, we find 



/ 



dESfi = e 



C a 



C ift e 3 v iF 



5Ui 



C 



ev iF 

(SUi -SUj). (A6) 



In the following, we write the distribution function as 



Here, we introduced Xi — ^lr(^l — ^r), the total ca- 
pacitance of cavity i, = J2 r Ci r + Ci g , as well as 
its electrochemical capacitance Ci^ = e 2 Vi F Ci^/(e 2 Vi F + 

Using Eq. (A6) to eliminate Sfi from the energy- 
integrated kinetic equation (Al), we obtain a set of cou- 
pled, nonlinear Langevin equations that determine SUi, 



(C iE + C)SUi - CSU-i = - G ii: 



c 



1 9 

e 2 v iF 



SU i + G' iS ^(SU i ) 2 + SI iS 



•Ci 



%\1 



(SUi - SUj) + G' t 



c 



[(SUi) 2 - SUiSUi 



(A7) 
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Appendix B: Diffusion coefficients 

The diffusion in Eq. (A7) is characterized by the diffusion coefficients defined as 37 

(SI ir (t)6I ir (0)) = ^ J dET ir [f ir (l - f ir ) + fi(l - fi) + (1 - T ir )(f ir - hf] 5{t) = D ir S(t). 

I 



(Bl) 



Importantly, the diffusion coefficients D ir depend them- 
selves on SUi through the energy-dependence of the 
transmissions T ir . This leads to a certain ambiguity when 
converting the Langevin equation into a Fokker-Planck 
equation, see below. Evaluating the above integral and 
expanding the diffusion coefficient to linear order in the 
applied voltage, we obtain 



D ir = ±k B ® i {G ir -G' ir 5U i ). 



(B2) 



Appendix C: Fokker-Planck equation 

Given a nonlinear Langevin equation of the form 

Xi = /i(x) + gij(*)Vj(i) 



(CI) 



where x = (#i,£2> • • •) and rjj(t) is a noise source satis- 
fying (rjj(t)) = and (Vi(t)Vj(t f )} = Sij5(t — t r ), one can 
show 38 that it is equivalent to a Fokker-Planck equation 
of the form 



dP 

~dt 



_d_ 



fi 



dgu 

— 9ki 
dx k 



l_d_ 

2~dx, 



(guQjiP) 



(C2) 

where Einstein's sum convention is implied. The param- 
eter a takes the values in the Ito prescription, 1/2 in 
the Stratonovich prescription and a = 1 in the Klimon- 
tovich prescription. 35 In our analysis, it turns out that 
only the Klimontovich prescription gives vanishing cur- 
rents in global thermal equilibrium. In our problem, we 
have gij = Sijgi such that the Fokker-Planck equation 
simplifies to 



dP 

~dt 



_d_ 

dx. 



dgi 



ld_ 



(C3) 



Multiplying the Fokker-Planck equation with Xk and 
XkXi, respectively, and integrating over all variables 
we obtain the following equations for the expectation val- 



ues 



d , , , n . _ / dg k 



d_ 
Jt 



dt 



\ dx k 9k 
(xkxi) = (xifk+Xkfi) +a{ri 
hi (gl)- 



(C4) 



dx. 



~9k 



dgi 
x k ——gi 
oxi 



(C5) 



To make a closer connection to the discussion above, we 
introduce gi = \J 2Di with the diffusion constants Di and 
obtain 



d_ 
dt 



(xk) = (fk) +ol • 



dx k 



(C6) 



dD k 



dt 



(x k xi) = (xif k + x k fi) +a(xi h x k 



dx k 



ODi 

dxi 



2S ki (D k 



(C7) 



By comparing the equation for (x k ) with the original 
Langevin equation, we furthermore obtain for the expec- 
tation value of the random currents 



(SI ir ) = ( y2D ir r) ir (t) ) = a 



ODjr 

dxj 



(C8) 



For the calculation of heat current, we also need cor- 
relators of the form (xigij(x)rjj(t)). In order to obtain 
them, we multiply the Langevin equation (CI) by x^, 



1 d 

2dt 



xf = Xifi(x) + x i g ij (x)r] j (t). 



(C9) 



Taking expectation values and equating with Eq. (C7), 
we find 

dD ir \ 



{Xidl ir ) = {Xi\j 2D ir 7] ir {t)) ~{D ir ) + Ot (xi~ ^ 



(CIO) 



Appendix D: Charge and heat currents 



by 



The charge current between lead r and dot i is given 



J dE(T° r - eT^SUiMr - /) + 5I ir (Dl) 



which can be rewritten as 
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, ir^ir j^j- T ^ ^ ^ ir^ir 

J-ir 



-(Vir-Vir) 



Gn 



(Vir-Vir^Ui-liGir-GldUi) 



Gn 
Gi, 



C 



e 2 Vi¥ 



SUi 



C 



e 2 Vi¥ 



SU- 4 



]} 



-5I ir . (D2) 



The expectation value of the fluctuating current part is obtained from Eq. (C8). 

The heat current between lead r and cavity i is given by Ji r = (Ui — Vi r )Ii r . Without an applied bias voltage, we 
have Jir = 5UiL r and, hence, 



{Jii 



Gi r 



c 



e 2 Vi¥ 



(mr) 



c 



(SUiSUi) {0) 



(6Ui5I ir )(°\ 



(D3) 



where the potential-current correlator can be obtained using Eq. (CIO). The superscript (0) on the expectation values 
indicates that they have to be evaluated to zeroth order in the applied bias voltage. 

The heat current up to linear order is again given by Ji r = 5UiIi r as both the nonfluctuating part of Ui as well as 
the expectation value of Ii r are of first order in the bias voltage. We find 



{Jir) 



G tr G iP{§Ui)(0){Vir _ v ^ _ G^ {Vir _ ViM5Ui)2){0) 



GiY> 

— Gi r 



GjY 

Gi,, 



C 
e 2 v i¥ 



GiY> 

(m) 2 ) {1) 



c 



e 2 u iF 



(5Ui5Ui)W 



(SUi6I ir )M. 



(D4) 



Appendix E: Onsager relations 

According to Onsager 39,40 the linear response coeffi- 
cients of charge and heat currents as a response to bias 
voltage and thermal gradients are related to each other. 
To verify the Onsager relation for our system, we expand 
the charge current through cavity 1 and the heat current 
between the two cavities to linear order in the bias AV 
applied to cavity 1 and the temperature difference A6 



between the reservoirs of the two cavities, 

(Iil) = GAV + £AO, (El) 
( J 2 s) = MAV + NAB. (E2) 

In the main text we already found that L = &bAilr/t#c. 
Evaluating similarly the heat current in response to an 
applied bias voltage, we find M = — &bOAilr/t#c in 
agreement with the Onsager relation L = — M/0. 
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